Identification of the Orbital Pairing Symmetry in UPt 3 
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This paper summarizes the results of a comprehensive analysis of the thermodynamic and trans- 
port data for the superconducting phases of UPt3. Calculations of the transverse sound attenuation 
as a function of temperature, frequency, polarization, and disorder are presented for the leading 
models of the superconducting order parameter. Measurements of the specific heat, thermal con- 
ductivity, and transverse sound attenuation place strong constraints on the orbital symmetry of 
the superconducting order parameter. We show that the superconducting A and B phases are in 
excellent agreement with pairing states belonging to the odd-parity E2u orbital representation. 



PACS numbers: 74.25.Ld, 74.25. Fy, 71.27.+a 



I. INTRODUCTION 

Unconventional superconductivity, the electronic ana- 
log of superfluidity in 3 He, was discovered in the heavy- 
fcrmian metals UBei3 and UPt3 more than a decade 
ago.lira As in liquid 3 He the observation of multiple 
superconducting phases was the direct evidence foe. a 
multi-component superconducting order parameter.ETTJ 
The phases of UPt3 have since become a paradigm for 
unconventional superconductivity. However, unlike the 
case of 3 He the identification of the orbital and spin sym- 
metry of the order parameter has been a more difficult 
task. Heavy fermion metals are more complex materials 
in which disorder, magnetism, spin-orbit coupling and 
anisotropy must be factored into any realistic theory of 
superconductivity in these systems (c.f. Refs. |?],||). 

In this paper we present new theoretical results and 
analysis of the transport properties of the leading mod- 
els for the superconducting phases of UPt3 . These mod- 
els yield qualitatively different predictions for the trans- 
port properties in the superconducting phases. We cal- 
culate the ultrasonic attenuation for the A and B phases 
and discuss its sensitivity to order parameter symme- 
try, polarization direction and disorder. From .our anal- 
ysis of experimental data for the heat capacity ,EI thermal 
conductivityE30 and transverse sound attenuationEj we 
determine the topology of the excitation gap on the Fermi 
surface and conclude that the orbital symmetry of the or- 
der parameter in the A and B phases of UPt3 belongs to 
an odd-parity E2 U representation. 



II. PAIRING SYMMETRY 



ThCj-. -discoveries of multiple superconducting 
phasesaO'ESEJ of UPt3 led to several theoretical models 
for the superconducting phase diagram based on different 
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symmetry groups, or symmetry breaking scenarios Jl3~E3 
One class of models is based on a two-dimensional ( "E" ) 
representation of the hexagonal point group, D^h, with 
the multi-component superconducting order parameter 
coupled to a symmetry breaking Held (SBF). There are 
four E-representations for strong spin-orbit coupling: 
two E-reps for both even-parity (Ei 9 , E2 ff ) and odd- 
parity (Ei„, E2 U ) pairing. The E-rep models require a 
weak SBF that lowers the symmetry of the normal state, 
splits the superconducting transition and produces mul- 
tiple superconducting phases. cB The SBF is generally 
assumed to be the in-plane antiferrouiagnetic order pa- 
rameter that onsets at Tm ~ 5 A'jHJ however _pther 
explanations of the SBF have been suggested.EirEil The 
precise structure of the short-range AFM correlations, 
e.g. the spatial structure of domains, as well as the role 
of AFM as. a SBF for superconductivity is still an open 
questionEa One of the outcomes of the calculations sum- 
marized below is that simple model of equal-size, equally 
populated multi-domain structures for the SBF is in dis- 
agreement with the anisotropy of the sound attenuation. 

The models that have been most successful in explain- 
ing the properties of the superconducting phases of UPt3 
are based on the even-parity (E lg ) and the odd-parity 
(E 2 . u ) representations of the hexagonal point group. The 
Ei s representation is a realization of spin-singlet, d-wave 
pairing for a metal with a uniaxial symmetry, while the 
E211 model describes the hexagonal analog of spin-triplet, 
f-wave pairing. These pairing states have an orbital or- 
der parameter of the form A(p/) = 771^1 (p/) + 772^2 (P/), 
where 3^1,2 (P/) are the basis functions for the appropri- 
ate E-representation, and the amplitudes 77 = (771,772) 
transform as a two-component 'vector' under the same E- 
representation. Thus, in the E-representations the order 
parameter of the A phase is identified as r\ = (1,0), the 
B phase as rj = (1, i), and the C phase as rj = (0, 1) (see 
Fig. [T]) . These identifications then refer to the specific ba- 
sis functions, 3^i,2(p/), f° r a particular E-representation 
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given in Table M. The orbital order parameter differs sig- 
nificantly for the two models, particularly for the high 
temperature A phase (772 = 0). For Ei g pairing the A 
phase has the structure, A^ ~ PzPx, which has an equa- 
torial line node in the basal plane, as well as a longi- 
tudinal line node circumscribing the Fermi surface. For 
the E2« representation, ~ p z (j? x — p y ) also has an 
equatorial line node, but has two longitudinal line nodes 
oriented 90 degrees to one another. The low-temperature 
B phase of both models breaks time-reversal symmetry 
(with 772 ~ Hrji). As a result the longitudinal line nodes 
are closed by the growth of the second component of the 
order parameter; and for T — » 0, A# ~ Pzipx + iPy) for 
Ei g symmetry, while A B ~ Pz{Px + ip y ) 2 for the E 2u 
representation. Thus, the low energy excitation spec- 
tra for the B phase of the Ei 3 and E2„ models is de- 
scribed by an equatorial line of zero energy excitations 
(p z = 0) and pairs of point nodes of the excitation gap 
{Px — Py = 0) on the Fermi surface. There is a slight 
difference in the density of states from the point nodes 
because the gap varies linearly near the point nodes for 
the Ei s model, |A(p/)| ~ |t?|, but quadratically for the 
E 2u model, |A(p/)| ~ |t?| 2 . These slight differences are 
predicted to be ob^r-vable in the heat transport at ultra- 
low temperatures .Ea~E3 




T 

FIG. 1. The phase diagram of UPt3. The three supercon- 
ducting phases A, B, and C with amplitudes r\ — (771, 772) meet 
with the normal state (N) at the tetracritical point. For sim- 
plicity the additional two Meissner phases are not shown. 

All of the E-rep models are based on two-component 
orbital order parameters. However, they yield different 
predictions for the thermodynamic, magnetic and trans- 
port properties, including the H-T phase diagram. One 
important difference arises for the case of weak in-plane 
hexagonal anisotropy, as is-.reflected by the very small in- 
plane anisotropy of i? c2 .c3E3 Weak in-plane anisotropy 
leads to an apparent tetracritical point for all field ori- 
entations provided the_prder parameter belongs to an E2 
orbital representation.!] 

A key difference between the even-parity and odd- 
parity E-representations is the spin structure of the or- 
der parameter. Odd-parity representations are spin- 



triplet pairing states, and in the absence of spin-orbit 
coupling the dimensionality of these representations is 
three times larger than that of the corresponding spin- 
singlet E-representations. However, strong spin-orbit 
coupling in the Uranium-based heavy fermion metals re- 
duces the symmetry group by allowing only joint rota- 
tions of the spin and orbital degrees of freedom. The 
even and odd parity representations are still described 
by (pseudo) spkbsiaglet and spin-triplet order parame- 
ters of the forrr£3~E3 

A Q/3 (p/) = A(p/) (ia y ) a0 (singlet) , (1) 
Aa/3(P/) = A(p/) • {i(T(T y ) al 3 (triplet) . (2) 

The triplet representations transform only under joint 
spin and orbital rotations of the discrete point group for 
the normal state, i.e., A(p/) — ► 7£A(7?. -1 p/), where the 
rotation 1Z G [D6h] tpi „_ OIbit . The full symmetry group of 
the normal state is Q = [-D6/i]st>m-orbit x T x U(l) with 
[-C>6h] S pin-orbit representing the hexagonal point group with 
inversion, T is the time-inversion operation and U(l) is 
the group of gauge transformations. In the limit of no 
spin-orbit coupling A(p/) transforms as a spin vector 
under the vector representation of the full spin-rotation 
group, and separately as a representation of the point 
group with respect to the orbital momentum, pf, i.e., 
A(p/) -> ftspmA^-j^p/), where TZ spin € SU(2) spin and 
7^-orbit G [Dehlaibit- In the absence of spin-orbit coupling 
the enlarged symmetry group for the normal state is 
G = SU{2) apia x [D 6h ] OTbit x T x U(l). 

There are two special classes of spin-triplet order pa- 
rameters that are frequently discussed as candidates for 
the phases of UPt3 . The first class are states in which the 
spin-triplet order parameter factorizes into a single spin- 
vector and an orbital amplitude, i.e., A(py) = dA(pj) 
where d is a real unit vector and A(p/) is an odd-parity 
orbital function. The vector d defines the axis along 
which the pairs have zero spin projection, e.g., if d||z 
then A TT = A lL = and A u = A iT = A(p). If we 
choose the quantization axis to be perpendicular to d, 
i.e., d _L z, then the same pairing state is described as 
equal spin pairing in an "easy-plane" , i.e., the pairs form 
triplet states with amplitudes A<- = A^ = A(p) and 
A^ = 0. In the second class, d is complex and the 
spin components of the order parameter spontaneously 
break time-reversal symmetry. In the general case A is 
complex, with A x A* 0, and varies over the Fermi 
surface. These states are called 'non-unitary' because 
the square of the spin-matrix representation of the or- 
der parameter is no longer proportional to the unit spin 
matrix, [A+A]^ = |A| 2 5 a/3 + i [A x A* • cr] a[j . As a 
consequence the spin degeneracy of the excitation spec- 
trum is lifted and the quasiparticle energy depends on 
the local pair spin at p/: S pair (p/) ~ iA(p/) x A(p/)*. 

Whether or not spin-orbit coupling is weak or strong on 
the energy scale of ksT c has important implications for 
both the orbital and spin components of thp-order param- 
eter that are allowed by symmetry. BlountH and Volovik 
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and GorkovEII showed that line nodes are not required for 
odd-parity states when spin-orbit coupling is relevant. 
However, line nodes in the a6-plane of the Fermi surface 
are allowed, and required for some representations, if the 
normal-state spin-orbit interactions lock d along the c 
axis of the crystal, i.e., d||c. Precisely this orientation of 
d was predictedH for UPt3 based on anisotropic param- 
agnetic limiting.Ej This effect arises from the competition 
between the condensation energy and the Zeeman energy. 
For d locked along the c axis of the lattice the Zeeman 
energy is pair-breaking for H||c, giving rise to paramag- 
netic limiting. However, for H 1 c the Zeeman energy, 
•^zeeman ~ (d • H) 2 , is minimum (vanishes); as a result 
there is no paramagnetic limit for this field orientation. 
The anisotropic paramagnetic limiting of H C 2 is sensitive 
to the spin structure of the order parameter , .but insen- 
sitive to the orbital pairing symmetrySoQO The odd- 
parity Fi2 U representation with strong spin-orbit locking 
of d||c quantitatively accounts for the anisotropy of the 
paramagnetic limit of H C 2 observed at low temperatures. 

The spin-singlet Ei s model appears to be incompat- 
ible with both the tetracritical point for H JL c and 
the anisotropic .paramagnetic limiting of H C 2- However, 
Park and JoyntE3 argue that there is enough freedom in 
the Ei 9 model to account for the existing experimental 
data on H C 2 . Both E-rep models have recently been chal- 
lenged by observations of ^nearly temperature indepen- 
dent Knight shift for H| |c,E3 which is interpreted in terms 
of non-unitary_spin-triplet pairing with weak, or no spin- 
orbit couplingH The authors of Ref. [37] assume that the 
Knight shift measures the bulk spin susceptibility. If, for 
simplicity, we ignore the anisotropy of the normal-state 
susceptibility, then for a given orientation of d the spin 
susceptibility is given by Xij = Xn (% - ckdj) + Xo didj, 
where \n is the normal state spin susceptibility and 
Xo(T) is the spin susceptibility for H||d, which is sup- 
pressed by pair-breaking and vanishes for T — * in the 
clean limit. For strong spin-orbit coupling with d locked 
along c we expect a suppression of the Knight shift for 
T < T c for fields H||c, but no suppression for H 1 c. 



However, in the limit of no spin-orbit coupling the Zee- 
man energy is minimized by rotation of d perpendicu- 
lar to the field. This implies that the Knight shift will 
be temperature independent and given by the normal- 
state shift for all field orientations. The NMR measure- 
ments of the Knight shifttZl appear to be in conflict with 
anisotropic paramagnetic limiting of H C 2- The paramag- 
netic limit observed for h]} 2 is a robust, thermodynamic 
property of bulk single crystals of UPt3, and a consis- 
tent interpretation of the NMR results for the Knight 
shift data must accommodate anisotropic paramagnetic 
limiting. This cannot be accomplished with a model of 
spin-triplet pairing without strong spin-orbit coupling. 

Our analysis presented below for the heat capacity and 
transport measurements is independent of the interpreta- 
tion of the Knight shift measurements. We show that the 
heat capacity, low temperature thermal conductivity and 
transverse sound attenuation data, in addition to the H- 
T phase diagram, are in quantitative agreement only for 
the odd-parity F,2 U representation (Table |) , independent 
of the orientation of d. In order to demonstrate this fact 
we present calculations for other models that have been 
proposed to account for the phase diagram. Thus, in ad- 
dition to the 2D E-representations we also examine the 
transport properties of the order parameter models be- 
longing to mixed r epres entations of the -D^r-point group, 
i.e., the AB modelsMnla and the AE model.E^I These mod- 
els were proposed as alternatives to the E-representations 
to explain the Ginzburg-Landau region of the H-T phase 
diagram. The most promising candidate of the AB model 
is the odd parity, spin-triplet model with mixed A2 U ©Bi u 
symmetry. The orbital order parameter for the A phase 
has the form, A^(p^) ~ p z rm(p x + ip y ) 6 , exhibiting an 
equatorial line node and six longitudinal line nodes. We 
also analyze the transport properties of the even-parity, 
spin-singlet Ai s ©Ei s model with an A phase of the form, 
Ayt(p/) ~ (2p 2 — p x — Py), which has a pair of 'tropical' 
line nodes located off the equatorial plane. For a more de- 
tailed description of the order parameter for these models 
see Refs. 



TABLE I. Polynomial functions representing the symmetry of the low-temperature B phases of several pairing models. The 
first three entries are based on the symmetry group [Ds^spm-orbu x T x (7(1). The third entry is representative of the class 
of AB models, and the last entry belongs to mixed symmetry representations resulting from the crystal-field splitting of the 
enlarged symmetry group, SO(3) 3pin - orbi t x T x 17(1). 

r yr point nodes line nodes cross nodes 

El g Pz{p X +iPy) tf = 0,7T 0=f 

E 2u Pz(Px+ip y ) 2 tf = 0,7r tf=f 

A2 U © Biu Ap z lm(p x + ip y ) 6 ip n =n^, = 0, 7r A <£>n 

-MBIm^ + ip y ) 3 n = 0, ..,5 

Ai g ®Ei g A (2p 2 z — pi— p y ) 1 9 = cos- 1 ±i 

+ i Epyp z A = 0, n 
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III. TRANSPORT THEORY 

Electronic transport in the superconducting state is 
sensitive to the nodal structure of the oidetnr^acam- 
eter, A(p/). Recent theoretical analysescIOO'E!! of 
low-tempcratu^e-thermal conductivity data on supercon- 
ducting UPtaEl'EJ have eliminated most of the theoret- 
ical models proposed to explain the phase diagram of 
UPt3. The non-unitary, spin-triplet pairing states based 
on a one-dim^iisipnal (ID) orbital representation stud- 
ied so far,c3E3Ea as well as the two-component order 
parameter models obtained icom nearly degenerate one- 
dimensional representationsEj ('AB models'), are unable 
to describe, even qualitatively, the temperature depen- 
dence and anisotropy of the thermal conductivity at low 
temperatures. The only pairing models which can ac- 
count for the thermal conductivity data are the two- 
dimensional (2D) orbital representations, Ei 5 and E2 M , 
and the Ai g ©Ei g (AE) model. However, the AE model 
predicts a large afy-plane anisotropy, which has so far 
not been observed.c2l After it was shown that a non- 
unitary, spin-triplet state with a ID orbital basis function 
was incompatible with the thermal conductivity data,c2l 
Machida et alo modified their weak spin-orbit coupling 
model by adopting the 2D orbital representation E 2ll . 
However, the model of Ref. proposes a spin structure 
for the order parameter which is in conflict with the ob- 
served Pauli limiting for H||c, and it predicts a fourth 
superconducting phase which disagrees with the phase 
diagram. 

Transverse ultrasound is an even more powerful probe 
of the order parameter and excitation spectrum than 
the thermal conductivity.elrEj The attenuation of hy- 
drodynamic sound is determined by the electronic vis- 
cosity tensor, which is sensitive to the relative orienta- 
tion of the polarization, propagation direction and or- 
der parameter r 8 lr 9 l The broken symmetries of the pairing 
state give rise to additional anisotropy of the low-energy 
excitation spectrum that is specific to the pairing state; 
the selection rules for acoustic absorption reflect these 
broken symmetries. 

In the hydrodynamic limit, lot <C 1 and q£ -C 1, where 
£ = VfT is the quasiparticle mean- free-path and r is the 
transport collision time, the ultrasonic attenuation is de- 
termined by components of the viscosity tensor 



a(q,£,T) = (lo 2 / gc 3 s ) n i: j ikl (q,Lo) £ l q ] i k qi , 



(3) 



where g is the mass density, c s — lo / q is the speed of. 
sound mode with wavevector q and polarization e£! 
The hydrodynamic limit is achieved even for high-purity 
single-crystals of UPt3. For the experiments reported 
in Ref. |l2| with propagation q||a and polarization e||b 
the sound frequency \&-io/2ir ~ 165 MHz, the speed of 
sound is c s ~ 2.1 km/sj23 and the elastic meanJijee path 
is l a b — Vf.abT ~ 1.5km/s • 240 ps ~ 360nm,E3, yield- 
ing lot « 0.25 and q£ a b ~ 0.18 at T = 0. Similarly, for 
q||a and e||c the reported values are lu/2it ~ 228 MHz, 



c s ~ 1.4km/s, and £ c « y/27f£ ao , yielding lot k, 0.34 and 
q£ c « 0.61 at T = 0. The parameters q£ c and lot are a 
factor of two smaller near T c than they are at low temper- 
ature, since t(T c ) ~ r(0)/2. Nevertheless, attenuation 
measurements for e||c are near the borderline of the hy- 
drodynamic regime. Measurements above and below this 
cross-over regime would be desirable; both for checking 
the applicability of hydrodynamic results for the atten- 
uation for e||c, and such measurements might exhibit to 
new phenomena in the collisionless regime. 

The viscosity and sound attenuation are calculated 
from the response of the momentum stress tensor to an 
ionic displacement field A(q, lo) = A(q,to)e. For trans- 
verse modes (q • e — 0) the stress and v-isco^ity tensors 
are related in the hydrodynamic limit by 53 ' 



Ily (q, to) = lo q A(q, lo) r) ijt hi (q, to) £ k qi 



(4) 



At low temperatures the transfer of energy and mo- 
mentum between the ionic lattice and electronic excita- 
tions is dominated by the scattering of quasiparticles off 
impurities or defects. The theory of momentum trans- 
port by quasiparticle scattering is formulated in terms of 
nonequilibrium Green's functions for electronic quasipar- 
ticles coupled to the acoustic modes of the lattice. The 
momentum stress tensor is 



n«(q,w) =N f J-^- Jd Pf [v / ] i tp/]i<5fl^ r (p/,q;e,w), 



(5) 



where v/ is the Fermi velocity, p/ is the Fermi momen- 
tum, Nf the density of states at the Fermi surface, and 
6g is the nonequilibrium quasiparticle Green's func- 
tion, integrated with respect to the quasiparticle energy, 
£ p ~ Vf(p — pf); 5g K includes both the changes in the 
distribution of occupied states and the dynamics of the 
spectrum of low-energy excitations (see Appendix) . 



IV. TRANSPORT EQUATIONS 

We use Keldysh's formulation of the nonequilibrium 
response theory and calculate the transport properties 
and related Green's functions in the quasiclassical limit, 
which is easily achieved in UPt3 for excitation energies, 
(fceT, Tilo) <C Ef w 1 meV, and wavelengths long com- 
pared to the Fermi wavelength, i.e., Tiq <ti pf. The cen- 
tral equation for the nonequilibrium Green's function in 
the quasiclassical limit is a transport equation. For small 
deviations from equilibrium the transport equation may 
be linearized in the deviations of the Green's function 
from its local equilibrium form. Our analysis and nota- 
tion follows that of Ref. |d], which provides a detailed 
discussion of the quasiclassical linear response theory, in- 
cluding a complete solution to the linearized nonequi- 
librium transport equations. A summary of these equa- 
tions, applicable to momentum transport in Fermi-liquid 
superconductors, is given in the Appendix. 
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The transport equation for the Green's function in- 
cludes the acceleration of electronic quasiparticles by the 
acoustic field and collision terms which transfer momen- 
tum between the lattice and the electrons. The stress ten- 
sor, and therefore the electronic viscosity which damps 
the acoustic wave, is calculated from the solution of the 
transport equation for the nonequilibrium Green's func- 
tion, Sg , which is driven by the coupling of quasiparti- 
cles to the ionic displacement field, i.e., an externally_im=. 
posed sound field, cr ext (p/,q;u) = i(yf -q)(p/ • A)i.E3o 
Below we report new results for the electronic shear vis- 
cosity for the order parameter models of UPt3, and new 
calculations of the anisotropy and temperature depen- 
dence of the attenuation which we use to interpret the 
experimental data for UPt3. 

In the limit of u> — ► and for resonant scattering the 
viscosity tensor simplifies to 



N 



f 



1 



d Pf 



&TT 3 k B T 

M?[P/]' run* 



J desech 2 (e/2k B T) x 



ReC R 



[9o So ~ Jo Jo + 11 \ 



(6) 



where C R = -^\A(p f )\i - (e~«) 2 , g R = i R /C R , f R = 
-A(p f )/C R , and e R = e - \Tv[T 3 a R p ] is the impurity- 
renormalized energy. In the case of triplet pairing with 
a unitary order parameter the only change in Eq. (|J) is 
replacement of A -> A and f R f R * -> • f fl *. See the 
Appendix for details on the notation. 

For a normal metal with a spherical Fermi surface 
we obtain Pippard's result for the viscosity rjij ij — 

_2_„,2 „2 
15 



jpjNfT for i 7^ ja Below T c the sound attenua- 
tion drops; for a conventional superconductor, in the 
limit C 1, rthe attenuation decreases exponentially 
for k B T -C A.E3 But for an unconventional supercon- 
ductor in which the order parameter vanishes at points 
or lines on the Fermi surface, the attenuation decreases 
with temperature as a power law reflecting the spec- 
trum of low-ener gy , exci tations near the nodes of the or- 
der parameter e3'c3E3£3 Impurity scattering modifies this 
spectrum near the nodes, and at low energies a new en- 
ergy scale, 7 <C Aq, appears, which is roughly the 'band- 
width' of low-energy quasiparticles bound to the impu- 
rity distribution by Andreev scatteringMil The bandwidth 
also appears as an impurity-renormalized quasiparticle 
width at zero energy, e R (Q) — ry. This new low energy 
scale defines a cross-over from the power law behavior as- 
sociated with scattering of continuum quasiparticles for 
"f/k B < T <^iT c , to a temperature ipdependent attenua- 
tion in the limit k B T < 7 HBEK3 

The bandwidth of the impurity-induced Andreev lev- 
els is determined by the self-consistency equation for the 
quasiparticle self-energy, 



t = r u ■ 



( 7 [|A( P/ )| 2 +7 2 ] 



cot 2 «5 + (7[|A( P/ )|2+72]-2): 



(7) 



where (...) is an average over the Fermi surface and 
r« = ^imp/ 7 !"-^/ is the scattering rate in the normal 
state for resonant impurities and 5o is the scattering 
phase shift. For the high purity UPt3 crystals studied 
in Ref. [l0|, i.e., low scattering rate, the crossover tem- 
perature is very low compared to T c ; analysis of the 
thermal conductivity provides a determination of both 
the scattering phase shift as well as the bandwidth of 
the impurity-induced Andreev states. The scattering 
centers are nearly resonant, i.e., 5q ~ t/2, giving a 
bandwidth, and cross-over temperature, of order 7 ~ 
k B T* w 0.2y/fiT u A fa 0.07fc B T c , where T u ~ 0.03k B T c , 
A — 2.Qk B T c and the slope of the excitation gap near 
the line node is fx = Aq 1 |dA/di?|^ =7r / 2 — 2. Thus, trans- 
port experiments on UPt3 have so far not investigated 
the ultra- low temperature region k B T <C 7 in any sys- 
tematic way. 



V. TRANSVERSE SOUND ATTENUATION 



The anisotropy and temperature dependence of the 
sound attenuation is sensitive to the polarization of the 
sound field and the symmetry of the order parameter. 
This was the basis of transverse sound attenuation ex- 
periments that provided early evidence for a line of nodal 
excitations in the basal plane.cll We examine the ab- 
plane anisotropy of the transverse sound attenuation. 
Our analysis covers the full temperature range below 
T c , and is particularly sensitive to the polarization and 
and anisotropy of the order parameter for both A and B 
phases of UPt3. To illustrate the sensitivity of the trans- 
verse sound polarization to the order parameter symme- 
try consider the theoretical models for the A phase of 
UPt3. The ab- plane anisotropics of the excitation gap, 
|A($ ^ 7r/2, <fi)\, for the A phase of four pairing models, 
E 2ii , Ei g , A 2u ® B lu and A lg E lg , are shown in Fig. §. 

The propagation and polarization vectors determine 
the angular dependence of the momentum transport by 
quasiparticles on the Fermi surface; the matrix element 
is proportional to [p/]f [ v /] 2 - This angular dependence is 
weighted by the angle-resolved density of states for mo- 
mentum transfer via impurity scattering, which depends 
on the anisotropy of the quasiparticle excitation spec- 
trum through |A(p/)| 2 . When both the propagation a,ud 
polarization vectors are in the basal plane (q||a, e||b)E3 
the matrix element is proportional to sin 2 2<f>, and is max- 
imum at angles of 7r/4 from these two axes, i.e., the mid- 
points between the polarization and propagation direc- 
tions. If these midpoint directions coincide with nodal 
directions (e.g. |A(p/)| ~ |cos20| for E 2u ) then the 
attenuation will be a maximum, while if the midpoint 
directions are along the antinodal directions then the at- 
tenuation is a minimum. This is illustrated in Fig. ^ 
where the polarization is directed along the a direction. 
The attenuation is largest when the polarization is along 
an antinode of the order parameter, and it is smallest 
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when the polarization is along a nodal direction. One 
can immediately see that we should expect to observe a 
rather different afo-plane angular dependence to the at- 
tenuation for the different order parameter models pro- 
posed for UPt3. We quantify these remarks below. 




FIG. 2. Sketch of the <j> dependence of the A phase pairing 
states at i? 7^ 7r/2 for (a) E 2u , (b) Ei g , (c) A 2u © B lu , and 
(d) Ai g (B Ei g models. The crossed lines represent the areas 
of maximal absorption of a sound probe with ab symmetry 
(q||a, e||b) at the Fermi surface. The dark shaded areas rep- 
resent the distribution of quasiparticle excitations. 

A. Results 

In order to make quantitative predictions for UPt3 
we use heat capacity and thermal conductivity measure- 
ments to fix the magnitude of the order parameter, the 
Fermi surface anisotropy, the nodal parameters and the 
scattering rate, all of which control the temperature de- 
pendence-arid anisotropy of the thermal conductivity be- 
low T c £3'c2l We used variational basis functions based on 
the polynomial functions in Table |. Symmetry dictates 
the geometry of the line and point nodes, as well as the 
topological indices for the point nodes, but not the slopes 
or curvature of the nodal regions of the excitation gap. 
Our variational procedure fits these parameters to the 
low-energy excitation spectrum from the thermal con- 
ductivity data of Ref. The order parameter is then 
determined self-consistently. For a detailed discussion of 
this analysis see Refs. |2^^. The anisotropy and tem- 
perature dependence of the transverse sound attenuation 
are then calculated with no additional parameters or ad- 
justments of the order parameter models shown in Fig. 

The transition temperature and splitting of the zero- 
field transition in UPt3 determine the scale and rel- 
ative magnitudes of the two order parameter compo- 
nents; they determine the dominant and subdominant 
instability temperatures. The instability temperature 
for the dominant channel is the transition temperature, 
T c i = T c+ , while the second instability temperature rep- 
resents the strength of the subdominant pairing channel, 
T C 2 oc uj c exp(— I/V2), where V2 is the pairing interac- 
tion for the, sub-dominant channel and lo c is the cut- 
off energy.tS The physical transition temperature, T c _, 
separating the A and B phases depends on the relative 
strength of the two paring channels, i.e., on T C 2/T C+ < 1. 
We solve the gap equation for the second transition and 
adjust T C 2 to the observed splitting Wrthc specific heat, 
(T c+ - T c -)/T c+ fa 60mK/495mKiiy All other mate- 



rial parameters are taken from out previous analysis of 
the thermal conductivity of UPt3.c9 We show the quality 
of the theoretical fits to the thermal conductivity, k, and 
heat capacity for the E2 U model in Fig. pi Similar fits can 
be obtained for the Ei g and AE models. However, only 
the in-plane thermal conductivity data can be accounted 
for by the AB model. 




FIG. 3. nfiit of the E2 U model to the thermal conductivity 
data (top)E3 and the specific heat data (bottom) B 

1. Identification of the A-phase 

In Fig. [| we show the attenuation data reported by 
Ellman et aZ.EHI for transverse sound propagation in the 
a6-plane with q||a and polarizations both in- and out- 
of-plane, e||b and e||c. These measurements were made 
on the same batch of crystals of UPt3 as the heat ca- 
pacity and thermal conductivity measurements shown in 
Fig. 2. In addition to the anisotropy associated with 
the polarization, the data show a pronounced change in 
the anisotropy and temperature dependence at the A-B 
phase transition. 

The enhanced absorption in the A phase compared to 
that of the B phase at the same temperature results from 
excess quasiparticles that scatter off the impurity distri- 
bution due to additional nodes of the A phase order pa- 
rameter compared to that of the B phase. Below T c _ the 
subdominant order parameter closes the additional nodes 
of the A phase order parameter. Thus, the sound absorp- 
tion drops faster in the B phase than it would in the A 
phase. The experimental results for the sound attenua- 
tion, including the anomaly at T c _ , are in excellent agree- 



G 



ment with theoretical calculations for the E2„ model of 
the A and B phases, with an A phase given by a rj = (1,0) 
state, but not a (0, 1) state; and a B phase at low temper 
ature which is (approximately) the r/ — (l,±i) state.EHl 
The comparison between theory and experiment is un- 
satisfactory for all other models. Even including the fre- 
quency dependence of the viscosity, ffa',^ (q, cj), does not 
change this result. We find only minor corrections to the 
sound attenuation for uj <C A (see Fig. ||). The (1,0) 
state of the A phase, as determined by sound attenuation, 
also agrees with the order parameter orientation obtained 
from the observed six- fold oscillations of H C 2 (<fr) , and the 
change in sign ofj-these oscillations when crossing the A-C 
phase boundaryEil 




FIG. 4. Comparison between the measured and calculated 
transverse sound attenuation of UPt3. The theoretical calcu- 
lation is for the E2« model and the phenomenological model 
previously used for the total elastic and inelastic scattering 
rate, F(T) = 0.03fc s T c+ (l + T 2 /T C 2 + ), obtained from the thep, 
mal conductivity analysis. The data are from Ellman et al.fj 
which are corrected to vanish at T = 0. 



In Fig. g we show the sound attenuation for ab and ac 
polarizations and different pairing models. Note in par- 
ticular that none of the order parameter models shown in 
Fig. |^ can account for both anomalies in a a b and a ac at 
T c _ . The main result of this work is that the experimen- 
tal data for the heat capacity, thermal conductivity, and 
sound attenuation, as well as the H-T phase diagram, 
are explained only by an order parameter with an orbital 
E211 pairing symmetry. We emphasize that there are no 
adjustable parameters in the calculation of the sound at- 
tenuation; all parameters of the model were previously 
determined by fitting the theoretical model parameters 
to the heat capacity and thermal conduct ivityHa 




0.0 0.5 0.0 0.5 1.0 



T/T T/T 

c+ c+ 

FIG. 5. Attenuation of ab and ac transverse sound for the 
various order parameter models. 



2. Domain Structure 



Neutron, diffraction studies in pure UPt3 under 
pressure,tJ and rat ambient pressure in Pd doped sam- 
ples (Pt «-» Pd) J23 demonstrated that the splitting of T c 
correlates with the basal plane AFM, suggesting a SBF 
coupling between the superconducting and AFM order. 
It has been argued, based on neutron diffraction studies 
as a function of magnetic field, that the AFM order is de- 
scribed either by a distribution of three equally populated 
domains, with Q vectop. oriented 120° to one another, or 
a triple-Q structure.Eji£j If a domain structure is present 
then, in the AFM model for the SBF, the superconduct- 
ing order parameter describing the A phase in the E-rep 
models may also form a domain structure. Such a do- 
main structure leads to a weakening of the anomaly for 
a a b at T" c _, but leaves the anomaly in a ac {T c ) virtually 
unchanged. 

In Fig. ^ we show the effect of multi-domain averaging 
on the anisotropy of the transverse sound attenuation. 
These calculations show that an AFM domain structure 
does not destroy the anomaly in the attenuation at T c _ . 
However, averaging over domains suppresses the charac- 
teristic enhancement of a a b for the E2 U order parameter 
coupled to a dominant or single domain of the SBF. In 
particular, the domain-averaged attenuation, (a a t>), for 
the E2« pairing state drops roughly twice as fast as the 
measured attenuation in the A-phase. If the AFM order 
parameter is the SBF for the superconducting phases, 
then our calculations are in agreement with transport 
and heat capacity measurements only for an E2 U order 
parameter coupled to a dominant domain or a triple-Q 
structure for the antiferromagnetic SBF.Eil 
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T/T 

c+ 

FIG. 6. Configuration averaged attenuation assuming three 
equally populated mono-domains for an Ei s and Ei2 U pairing 
state. For a comparison the results for the E 2u state of Fig. 
are re-plotted. 



3. Scattering phase shifts 

It has been pointed out in several studies that the 
temperature dependence of the transport coefficients in 
UPt3 is qualitatively consistent with strong scattering in 
the unitarito-Jimit. Our analysis of the thermal conduc- 
tivity dataJ13 and the attenuation dataEJ confirm that 
the scattering phase shift is near the resonant limit; 
from the .analysis of the thermal conductivity we obtain 
So > 80°,t3 while the transverse sound attenuation data 
implies a scattering phase shift Sq > 60°, as shown in 
Fig. |. 




FIG. 7. Sensitivity of the transverse sound attenuation to 
the scattering phase shift So = 30°, 60°, 80°, 90° for the E 2u 
pairing model and the same parameters as shown in Fig. 



4- Ultra-low temperature region 



The temperature dependence of the transverse sound 
attenuation is predicted to change qualitatively below the 
cross-over temperature, fc^T < 7. The zero-temperature 
limit is finite, reflecting the finite density of states at 
the Fermi level, and the leading temperature-dependent 
corrections are of the Sommerfeld type, 0[(T/j) 2 ] for 
ksT < 7. The limiting attenuation in the ultra-low tem- 
perature limit, ksT <C 7, is obtained from the viscous 
stress tensor for T — ► and uj — ► 0. In this limit the 
states contributing to the stress tensor in Eq. (^) are 
confined to the impurity-induced Andreev band of or- 
der 7 3> max(T, u>). The spectrum and self-energy are 
weakly energy dependent on this scale, thus, we can eval- 
uate the slowly varying parts of the integrand in Eq. (|^) 
at zero energy to obtain 



Vijy 



N f dp f [v/]f[ P/ L 2 



7 



|A( P/ )| 2 +7 2 P 



(8) 



dected the Sommerfeld corrections of 



where we have i 
order 0[(T/~/) 2 ]B 

Whether or not the zero-temperature limit is univer- 
sal, i.e., independent of the density and scattering cross- 
section for the impurities, depends on the polarization 
of the acoustic wave and the symmetry of the ground 
state. The damping is universal only for transverse 
waves propagating along the high symmetry directions 
in the a6-plane, i.e., for polarization e||b propagating 
along q||a, or vice versa. For these polarizations the 
ground state of all the models we discuss, except the 
AE-model, possess a universal limit for the attenuation. 
Table || summarizes the results for the Ei 9 and E2u 
models. For the ground state of either E-representation 
we find ?7ab,ab — VfPfNf/(8ixAo), where fx w 2 is the 
slope of the excitation gap near the line node in the 
a6-plane. For any other polarization the relevant vis- 
cosity is non-universal for u>,T -C 7. For example, 
for the E2 U ground state we obtain a limiting value for 
a *c(0) ~ ^A^ aab (°) <c "ab(7 1 c+ ), where [i 2 ~ 4 is the 
parameter defining the curvature of the excitation gap 
near the quadratic point node along the c axis.E3 

Further experiments using transverse sound can 
be used to confirm the predictions for the a6-plane 
anisotropy of the order parameter in UPt3. The ideal 
method would be to propagate transverse sound along 
the c axis and measure the attenuation as a function 
of the azimuthal orientation of the polarization in the 
ab-plane. The qualitative predictions for the in-plane 
anisotropy can be deduced from Fig. ||; a two-fold sym- 
metry of the anomaly at T c _ is expected for the Ei 3 
representation and a nearly isotropic attenuation above 
T c _ for the other models. In the B phase all models show 
a nearly isotropic attenuation except for the AE model, 
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which has a two-fold symmetry. A similar experiment 
was suggested by Moreno and ColemarO to map out the 
gap structure in the high-T c cuprates. 

TABLE II. Asymptotic low-temperature limits of the 
sound attenuation (lo, T — > 0) and thermal conductivity for 
the Ei s and E2u pairing states. The results are scaled in units 
of no = (n 2 /3)k 2 B TNfvfTA and a = (lo 2 /4gc 3 s )pf N/vfrA, 
where ta = ?i/2/iAo(0) is an effective transport time. 



Transport 
Coefficient 



E 



is 



^cc/^0 
Oab/QO 

Qac/flO 



2 M7 /(MiAo) 



8/j. Tu* 
1 + 2/i 2 Ao 



1 

1 

fi 2 A 



* In the strong scattering limit including vertex corrections. 

To summarize the main conclusions, we have shown 
that transverse ultrasound provides detailed information 
on the orbital pairing symmetry of both the supercon- 
ducting A and B phases of UPt3 ; the anisotropy and the 
anomalies in the temperature dependence of the attenu- 
ation for different polarizations is explained only by an 
order parameter. Measurements of ev a b and a &c are 
in excellent agreement with a (1,0) state for the A phase, 
corresponding to ap z (p 2 , — p y ) order parameter. Further 
measurements at lower temperatures, or as a function of 
impurity disorder, may also be used to test the predic- 
tion of a universal limit for the in-plane transverse sound 
attenuation. 
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APPENDIX A: 

In this appendix we summarize the relevant nonequi- 
librium transport equations for our calculations of the 
sound attenuation in unconventional superconductors. 
For a more extensive review of nonequilibrium transport 
theory in superconductors see Refs. f4l| , |68| . 

The nonequilibrium (Keldysh) Green's function, 



8g K = Sg H o $ - $ o Sg A + Sg a , 



(AI) 



contains the spectral response given by the retarded and 
advanced Green's functions (8g R ' A ), and the anomalous 



response given by 8g a , which in normal metals is essen- 
tially the nonequilibrium distribution function. The equi- 
librium distribution function is 4>o = tanh(e/2/csT), and 
we use the shorthand notation for the shifting product, 
$o o A = <f> (e - uj/2)A(e, uj) and A o $ = A(e, Lo)§ (e + 
lo/2). Pairing correlations and particle- hole coherence 
require a matrix structure for the particle-hole degree of 
freedom. The 'hat' over the Green's functions and self- 
energies indicates their 4x4 matrix structure in particle- 
hole and spin space. 

The transport equations for the Green's functions, lin- 
earized with respect to an external perturbation, are 



Sg 



R,A ;r,a 



~R.A 



o-oxt + 8 a 



R,A 



[(Text , $0] o 



g a oh J 
°9a~ 



Sa a og A 



9 R o 



[$0: 



g R o8a a 

<5ext] j 



(A2) 



(A3) 



where h RA — er^ — &q' A , ^3 is a Pauli matrix, and 
cr ox t is the external perturbation, e.g., the coupling of 
quasiparticles to a sound field (ionic displacement field). 
The transport equations determine the deviations of the 
Green's functions from their local equilibrium values, 
8q x = q x — q x , in terms of the external field and the 



corrections to the self-energies, 8a = a — a§ with 
X € {R, A, K}. The anomalous self-energy, 8a a , is de- 
fined similarly to the 8g a , 



8a 



K 



8a R o$ -$(|O Sa J 



(A4) 



The equilibrium Green's functions, g$ , and self-energies, 
a x , are inputs to the linearized transport equations. 
At low temperatures, and for long- wavelength, low- 
frequency sound, the damping of the acoustic wave is de- 
termined by the scattering of quasiparticles off impurities, 
and defects that are co-moving with the ionic latticed 
Impurity scattering enters the transport equations via 
the impurity-scattering self energies, 



X 



(p/,p/;e,t) , 



(A5) 



where the quasiparticle-impurity scattering t-matrices 
are given by 

i R ' A (p f ,p' f ; e, t) = u(p/,p'/) + N f (u(p f ,p'j) o 

g R > A (p'j; e, t) o f^p", p' /; e, t)) pl , , (A6) 

i K {p f ,p' f ;e, t) = N f (i R (p f ,p'}; e, t) o 

g K {p" } -e,t)oi A {p" f ,pi f ,e,t)) p , r (A7) 

The t-matrices are calculated self-consistently with the 
order parameter, A(p^), and with the impurity ver- 
tex, u(p/,pj), describing elastic coupling of quasipar- 
ticles to impurities. The calculations presented here as- 
sume isotropic, non-magnetic impurities with u = UqI 
for all (p/,p'|). This model is then described by two 







parameters, e.g., the s-wave scattering phase shift, So = 
tan -1 (ttN fUo) , and the density of impurities, n imp . An- 
other useful parametrization is in terms of the cross- 
section, g = (in/ft j) sin 2 So, and the normal-state scat- 
tering rate in the unitarity limit (So — tt/2), T„ = 
n^p/irNf, Note that the transport scattering rate in 
this model is given by 1/2to = T u sin 2 So, and the elastic 
mean free path is then £ cl = VfTQ. 

The other key term entering the transport equations is 
the off-diagonal pairing self energy, or order parameter. 
The general form for the pairing self energies is 

aR,a _ ( _ 0_ Aicr y + A ■ ia<T y \ , . > 

A - ^ Aia y + A ■ ia y a J ' (Ati) 

where the spin-singlet (A) and spin-triplet (A) order pa- 
rameters are given by the gap equations, 

A^(p /; t) = J ^-(VS(p f ,p' f )f K (pi f ;e,t)) py (A9) 

* R ' A ^) = j ^<Y t (p/,P^)-f^(p' / ;t)> p , , (Aio) 

where V s and V 4 are the pairing interactions in the sin- 
glet and triplet channels. The components A and A 
are related to A and A by fundamental symmetries (see 
Appendix C of Ref. ^); in equilibrium A = A* and 
A = A*. 

To complete the set of equations for the linear response 
equations we write the solutions for the equilibrium re- 
sponse functions in terms of renormalized quasiparticle 
energy and order parameter. For spin-singlet and unitary 
spin-triplet pairing the general solutions for the equilib- 
rium retarded and advanced Green's functions are 

~ R.A 

~R,A 6 T3 - A 

% = — 



(All) 



^R,A\2 _ ^R,A^2 



where the renormalized quasiparticle energy is 

e R ' A (p f ,e)=e-]Tr[f 3 a^(p f ,e)} . (A12) 



For resonant s-wave scattering 



-R,A, 



Pf 



(A13) 



For s-wave scattering the renormalization of the off- 
diagonal self-energy by impurity scattering vanishes, 



£*' (p/,e) = A(p/) ; 



(A14) 



for non- identity representations of the point group, i.e., 
excluding the A\ g representation. Finally, we note that 
the diagonal component of the nonequilibrium Green's 
function that determines the momentum stress tensor in 
Eq. (Q) is obtained from the matrix Green's function by 
8g K = \TrSg K . 
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